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1.  Summary 


We  present  a  theoretical  and  experimental  study  of  subpicosecond  kinetics  of  photo-excited 
carriers  in  bulk  gallium  nitride  (GaN).  The  theoretical  model  uses  a  simple  two-band 
approximation  with  the  conduction  band  nonparabolicity  taken  into  account.  The  scattering 
mechanisms  included  in  this  work  are  interactions  of  carriers  with  polar  optical  phonons  and 
carrier-carrier  Coulomb  interaction.  Interaction  of  the  semiconductor  with  the  exciting  laser 
pulse  is  treated  coherently,  while  the  scatterings  are  included  in  the  stochastic  Monte  Carlo 
solution  of  the  kinetic  equations.  The  theoretical  results  are  compared  with  the  experimental 
results  of  the  time-resolved  photoluminescence  induced  by  a  100-femtosecond  (fs)  pulse. 


2.  Introduction 


Ultrafast  phenomena  in  photo-excited  semiconductors  are  of  great  interest,  both  theoretically  and 
experimentally.  On  the  experimental  side,  ultrafast  laser  technology  improved  the  time  resolution 
of  measurements  of  carrier  dynamics.  On  the  time  scale  of  100  fs,  both  coherent  and  non¬ 
coherent  processes  contribute  to  the  observed  carrier  dynamics.  On  the  theoretical  side,  there  has 
been  much  progress  toward  explaining  the  experimental  results  related  to  dynamics  of  hot 
carriers  in  gallium  arsenide  (GaAs)  and  related  compounds  ( 1 ).  The  more  traditional  approach 
employs  the  semi-classical  Boltzmann  equation  that  describes  the  time  evolution  of  carrier 
density  distributions.  That  equation  can  be  studied  by  various  numerical  methods.  An  especially 
efficient  approach  to  finding  its  solutions  for  given  initial  conditions  is  the  Monte  Carlo 
simulation  (2).  Because  the  interband  polarization  is  not  included  in  the  Boltzmann  transport 
equations,  the  coherent  effects  due  to  the  interaction  of  the  external  light  field  with  the 
semiconductor,  and  the  excitonic  effects  due  to  electron-hole  interactions,  cannot  be  studied  in 
the  semiclassical  approximation.  For  the  purpose  of  modeling  various  ultrafast  phenomena,  the 
density  matrix  formalism  provides  a  consistent  description  (3).  This  approach  uses  perturbation 
theory  to  obtain  the  equations  of  motion  for  the  distribution  functions  of  hot  carriers  and  for  the 
interband  polarization,  and  takes  coherent  and  dephasing  processes  into  account  in  a  consistent 
way.  The  resulting  coupled  equations  are  referred  to  as  Semiconductor  Bloch  Equations  (3) 
(SBE),  and  they  replace  the  Boltzmann  transport  equation. 

In  good-quality  polar  semiconductors,  one  of  the  main  energy  and  momentum  relaxation 
mechanisms  in  carrier  kinetics  is  the  interaction  of  carriers  with  longitudinal  optical  (LO) 
phonons  via  the  Frohlich  interaction  (4).  In  GaAs,  the  relaxation  rate  of  electrons  at  slightly 
above  1  LO  phonon  energy  is  calculated  to  be  about  140  fs  (5).  GaN  has  larger  polarizability, 
and  the  electron-phonon  emission  rates  are  an  order  of  magnitude  higher  than  in  GaAs.  In 
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particular,  the  relaxation  rate  of  electrons  at  slightly  above  1  LO  phonon  energy  is  calculated  to 
be  about  12  fs  (6).  Therefore,  the  carrier  dynamics  in  GaN  can  be  substantially  different  from 
dynamics  in  GaAs.  This  feature  is  particularly  important  in  understanding  the  kinetics  of  the 
electron-hole  plasma  generated  by  sort  laser  pulses. 

The  carrier  kinetics  study  of  GaN  presented  here  was  motivated  by  the  current  interest  in 
application  of  nitride  semiconductors  in  ultraviolet  (UV)  optoelectronics,  and  high  temperature 
and  power  electronics.  The  study  can  be  extended  to  aluminum  gallium  nitride  (AlGaN) 
compounds.  The  UV-emitters  and  UV-detectors  in  the  wavelength  region  from  300-360  nm  and 
from  250-270  nm  will  serve  as  important  building  blocks  of  the  new  generation  of  biological 
agent  detection  systems.  GaN  is  potentially  important  in  many  applications,  such  as  power  and 
optoelectronics  devices,  because  it  has  a  number  of  favorable  features,  among  these  a  large  band 
gap,  high  saturation  velocity,  and  high  breakdown  electric  field.  A  better  understanding  of  hot 
carrier  dynamics  can  greatly  improve  the  knowledge  of  the  material  and  design  issues,  and  the 
performance  of  devices. 

In  this  work,  we  use  the  SBE  formalism  to  study  the  dynamics  of  carriers  in  bulk  GaN  under 
subpicosecond  laser  pulse  excitation,  and  we  compare  the  theoretical  results  with  the 
experimental  results  obtained  by  time-resolved  photoluminescence  (TRPL)  studies  (7).  The 
kinetic  problem  we  investigate  here  is  illustrated  schematically  in  figure  1 . 


Figure  1.  A  classical  light  pulse  whose  central  Fourier  component  has  frequency  v 

excites  electrons  and  holes  in  a  semiconductor  sample.  The  non-equilibrium 
distributions  of  electrons  in  the  conduction  band  and  heavy  holes  in  the  valence 
band  undergo  intraband  relaxation  and  subsequent  interband  recombination  with 
light  emission.  The  intraband  processes  included  in  this  work  are  carrier-phonon 
and  carrier-carrier  scatterings.  The  effects  of  the  trap  states  indicated  in  the  figure 
are  not  included. 

In  section  3,  we  compare  single-carrier  scattering  rates  due  to  various  interactions  and  compare 
their  thermal  averages  for  the  equilibrium  distributions.  We  find  that  with  the  material 
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parameters  established  for  the  bulk  GaN,  the  intraband  relaxation  is  dominated  by  the  carrier 
interactions  with  LO  phonons.  In  section  4,  we  establish  the  kinetic  equations — the  SBE  derived 
from  the  Heisenberg  equations  of  motion.  A  simple  relaxation  time  approximation  is  considered 
in  section  5,  where  the  effects  of  coherence  and  dephasing  in  the  interband  luminescence  are 
discussed.  The  theoretical  model  was  developed  in  conjunction  with  the  experimental  studies  of 
the  time-resolved  spectroscopy  of  UV  luminescence  outlined  in  section  6.  The  numerical 
solutions  of  the  SBE  were  performed  using  the  Monte  Carlo  method,  following  the  approach  of 
Kuhn  and  Rossi  to  account  for  the  effects  of  the  interband  polarization  (8,  9).  The  results 
corresponding  to  lower  excitation  power  are  presented  in  section  7,  where  they  are  compared  to 
the  available  experimental  results. 


3.  Intraband  Carrier  Scattering  Processes 


In  this  work,  we  use  an  isotropic  two-band  model  for  the  T  valley  of  conduction  electrons  with 
mass  me  =  0.2mo  and  heavy  holes  with  mass  mh  =  l.lmo,  where  mo  is  the  vacuum  electron  mass. 
The  nonparabolicity  of  the  conduction  band  of  GaN  will  be  taken  into  account  by  analytical 
approximation  similar  to  the  one  obtained  in  Kane’s  model  for  a  cubic  structure  (10) 

~~~  =  E(1  +  aE)  =  y(E)  ,  (1) 

2me 

where  a  =  (1  -  me/mo)  ( 1/Eg).  The  value  of  the  nonparabolicity  coefficient  a  obtained  for  the 
wurtzite  GaN  from  full  band  structure  calculations  (11)  turns  out  to  be  considerably  larger — 
about  0.37  eV  ’.As  an  approximation  here,  we  will  use  a  ~  1/Eg.  The  following  values  will  be 
assumed:  me  =  0.2mo,  mh  =  l.lmo,  and  Eg  =  3.39  eV. 

We  assume  dispersionless  LO  phonons  with  energy  hcoio  =  91.2  meV.  For  the  scattering  of 
carriers  by  LO  phonons  via  the  Frohlich  interaction,  the  matrix  element  in  a  sample  of  volume  V 
is 


M^(k,k  +  q) 


e  tmLO 

'  i  _ 1  y 

qss80V 

^£®  £J_ 

F(k,k  +  q) 


—12 

in  SI  units,  where  ss  =  5.35,  Soo  =  8.9,  and  the  vacuum  so  =  8.854x10  F/m.  The  overlap 
integrals  (4)  I(k.k')  are  given  by 


P  (k,k'j 


(l  +  °Et)'°(l  +  gE1-)'°+°(EtEt.)l°coseit. 
(l  +  2aEk  )  (l  +  2aEk, )] 


(2) 


(3) 


for  the  conduction  band,  and  by 
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Ih(k,k')  =  i(l  +  3cos2ekk,)1/2 


(4) 


for  the  valence  band,  where  9kk  is  the  angle  between  the  initial  and  final  momenta.  The 
scattering  rate  is  obtained  using  the  Fermi  Golden  Rule: 


3i 


r(k)  =  (27t//-l)vjTT|M(k,k')|2[(N1_k,+l)6(El-E1.-fltoLO)  +  Nk_l,5(Ek-El,  +  StoLO)]  (5) 


where  Nq  is  the  occupation  number  of  LO  phonons,  [exp(ha>Lo/kBT)  -  1]  ’. 


The  scattering  time  defined  as  h/T  is  shown  in  figure  2  as  a  function  of  carrier  energy  above  the 
band  edge  for  the  parabolic  case.  The  corresponding  momentum  and  energy  relaxation  times  ( 4 ) 
are  also  shown. 


Energy  (meV)  Energy  (meV)  Energy  (meV) 


Figure  2.  Left  panel:  scattering  times  for  electron  (solid  curve)  and  heavy  hole  (dashed  curve)  interacting 
with  polar  LO  phonons.  The  corresponding  momentum  (middle  panel)  and  energy  (right  panel) 
relaxation  times. 


In  the  presence  of  other  free  carriers,  the  carrier  interaction  with  polar  phonons  is  screened.  In 
the  static  limit,  the  screened  Frohlich  interaction  matrix  element  takes  form  (12) 


M^(k,k  +  q) 


e  hcoLO 

r  i  _ 

_  iv 

-1/2 

q 

8()V 

2  ,  2 
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where  the  static  reciprocal  screening  length  is  given  by  (3) 
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in  terms  of  electron  and  hole  distribution  functions  f6  and  fh.  The  effects  of  nonparabolicity  and 
screening  on  the  electron-LO  phonon  scattering  rates  are  shown  in  figures  3  and  4. 
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Figure  3.  Polar  optical  phonon  scattering  rate  for  non-parabolic  band. 
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II 


PO  Scattering  Rates  for  GaS  with  Screening 


Figure  4.  Polar  optical  phonon  scattering  rate  for  non-parabolic  band  with  screening. 

The  nonparabolicity  of  the  conduction  band  increases  the  emission  rate  of  LO  phonons  for 
electrons  with  energy  exceeding  0.4  V  above  the  band  edge.  This  is  what  prevents  the  runaway 
hot  electron  effect  in  GaN  in  strong  applied  field.  Both  the  nonparabolicity  and  screening  will  be 
taken  into  account  later  in  the  numerical  simulation  of  the  hot  carrier  kinetics. 
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The  carrier-carrier  scattering  proceeds  via  screened  Coulomb  interaction: 


1 


s0ssV  q  +  k" 


(8) 


with  the  screening  wave- vector  k  given  in  equation  7.  In  the  non-degenerate  case  of  carriers  with 
densities  ne  and  nh  at  temperature  T,  k  becomes  the  Debye  reciprocal  screening  length  defined  by 
k02  =  e2(ne+nh)/soSskBT,  while  in  the  strongly  degenerate  case  it  becomes  the  Thomas-Fermi 
wave-vector.  For  the  scattering  of  carriers  by  ionized  impurities,  we  assume  the  same  screened 
Coulomb  potential  as  in  equation  8,  the  Brooks-Herring  model  (4). 

The  criterion  for  the  applicability  of  classical  statistics  is  the  smallness  of  the  thermal  wavelength 
X  =  (27ih“/mekBT)  .  The  electron  gas  is  classical  if  X  n  «  1,  and  the  quantum  statistics  are  to  be 
used  if  X  n  >  1.  For  electrons  in  GaN  at  room  temperature,  we  obtain  kn  ~  0.9  if  n  =  10  cm  . 

1 8  —3 

Therefore,  the  electrons  become  degenerate  if  n  >10  cm  . 


Consider  next  the  scattering  of  carriers  by  acoustic  phonons.  For  the  deformation  potential 
interaction  the  matrix  element,  assuming  a  parabolic  band  and  linear  phonon  dispersion,  is  given 
by  (13) 


Mdp  (k,k  +  q)  =  D 


2puV 


(9) 


where  D  is  band  deformation  potential,  p  is  the  density,  and  u  is  the  sound  velocity.  For  GaN  p  = 

3 

6.15  g/cm  ,  u  =  6.56  Km/s,  and  D  =  8.3  eV  for  the  conduction  band.  For  the  unscreened 
piezoelectric  interaction 


MPE(k,  k  +  q) 


q2N)£s  L2puqV 


(10) 


where  e^j  is  the  third  rank  tensor  of  the  piezoelectric  constants,  2  is  a  unit  polarization  vector, 
and  summation  over  repeated  indices  is  implied. 

It  is  generally  the  case  in  polar  semiconductors  that  the  largest  contribution  to  the  electron 
scattering  rate  comes  from  the  electron  interaction  with  polar  optical  phonons  (3,  4).  In  order  to 
compare  the  scattering  rates  due  to  different  interactions  in  GaN,  let  us  consider  the  electron 
momentum  scattering  rates  (see  reference  4  for  the  definition)  and  then  evaluate  their  thermal 
averages  for  the  case  of  a  quasi-equilibrium  distribution  characterized  by  electron  temperature  Te 
for  the  non-degenerate  case  and  a  parabolic  band.  The  corresponding  momentum  relaxation 
times  are  shown  in  figure  5  as  functions  of  electron  temperature  Te/To  where  T0  =  300  K.  The 
electron  density  was  set  to  n  =  10  cm  ,  and  the  impurity  density  nimp  =10  cm  . 
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Figure  5.  Electron  momentum  relaxation  times  obtained 
from  thermally  averaged  momentum  scattering 
rates,  and  are  shown  as  functions  of  electron 
temperature  Te.  The  lattice  temperature  was  set  to 
T0  =  300  K,  electron  density  n  =  1018  cm"1. 


From  these  quasi-equilibrium  results,  we  expect  that  on  the  subpicosecond  scale  the  carrier 
kinetics  are  dominated  by  the  interactions  with  polar  optical  phonons.  In  the  carrier-carrier 
scattering  considered  previously,  we  assumed  static  screening.  When  the  dynamic  screening 
effects  are  taken  into  account,  one  finds  that  the  electron-hole  semiconductor  plasma  can  sustain 
collective  charge  density  waves,  or  the  plasma  oscillations.  An  electron  or  a  hole  whose  kinetic 
energy  exceeds  a  certain  threshold  value  can  lose  its  energy  to  the  collective  excitation.  The 
quantum  mechanical  picture  is  that  of  emission  of  plasmon  by  a  single  particle.  Alternatively, 
one  can  look  at  the  Fermi  quasi-particle  scattering  rate  evaluated  in  the  Random  Phase 
Approximation  (RPA)  and  obtain  a  similar  result  (14).  In  the  long  wavelength  limit,  the 
frequency  for  the  two  component  case  is  obtained  as 


co 


2 

0 


I 


e2n„ 


sns  m 

Us  o 


(11) 


2  2 

in  the  case  of  photo-excited  electrons  and  holes  ne  =  nh  =  n  and  coo'  =  ne  /mrsoSs,  where  we 
defined  the  reduced  mass  mr  =  (me_1  +  nth-1)-1.  For  GaN,  we  then  obtain 


/kon  =  30. 3x 


f  \ 

n 


1/2 


vnoy 


meV . 


(12) 


The  matrix  element  for  the  carrier-plasmon  interaction  is  given  by  (15) 


Ma-Pi(k,k  +  q): 


2*3 


t-fl 


2eoesma  v 


1/2 


2q  k  +  q2 

2®P  (q)q 


(13) 


7 


where  cop(q)  is  the  plasma  frequency.  Its  value  at  q  =  0  is  given  by  co0.  In  order  to  estimate  the 
electron-plasmon  scattering  rate,  we  consider  a  degenerate  electron  gas  characterized  by  its 
separate  Fermi  energy  Ep  =  (h  /2me)(37i'n)  ,  not  related  to  the  chemical  potential  of  the 
equilibrium  electron-hole  plasma.  The  plasmon  as  a  well-defined  collective  excitation  is 
restricted  to  q  <  qc,  where  qc  is  the  threshold  for  plasmon  decay  into  the  quasi-particle 
continuum,  obtained  from  the  zeroes  of  the  dielectric  function  in  RPA: 


qcM  ,  qc2ft 

me  2me 


(14) 


If  the  electron  energy,  Ek,  is  much  greater  than  fitoo,  the  rate  of  spontaneous  plasmon  emission 
can  be  approximated  (2)  as 


E(k) 


2  f  -M  A 

mee  ln  ftkqc 
4;is0es/rk  ^mero0y 


(15) 


which  can  be  rewritten  in  terms  of  electron  energy  E: 


f  _  V/4/c  V/2 


r  =  r, 


vnoy 


In 


^  E  ^ 


F 

vAoy 


E, 


xU2 


(r  V/2 


o  y 


F 

V^oy 


(16) 


1  Q  _ O 

where  1/T0  =  t()  ~  20.5  fs,  E0  =  hco0,  and  n0  =  10  cm  .  The  corresponding  scattering  time  is 
shown  in  figure  6  as  a  function  of  energy. 


Energy  (meV) 


Figure  6.  Electron-plasmon  scattering  time  for 

spontaneous  plasmon  emission  shown  as 
a  function  of  electron  energy  for  two  values 
of  density. 
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We  see  that  the  electron-plasmon  scattering  rate  is  by  a  factor  of  three  smaller  than  the  electron- 
LO  phonon  rate.  In  our  treatment  of  SBE,  the  dynamic  screening  effects  can  be  accounted  for  by 
some  of  the  terms  that  appear  in  the  second  order  in  carrier-carrier  interaction  (i).  We  did  not 
include  them  in  this  work,  but  they  should  be  included  in  the  future  work. 


4.  Derivation  of  Kinetic  Equations 


The  two-band  model  of  carriers  (conduction  band  electrons  and  valence  band  holes),  interacting 
with  phonons  and  subjected  to  a  semiclassical  electric  field,  can  be  defined  in  terms  of  creation 
and  annihilation  operators  for  electrons  and  holes  Ck+,  dk+,Ck,  dk  with  Fermi  anti-commutation 
relations,  and  for  phonons,  bk  and  bk+  with  Bose  commutation  relations  (3).  The  Hamiltonian  of 
the  system  is  given  by  H  =  Ho  +  Hc.Ph  +  Hc.c,  where  the  first  term  includes  non-interacting 
carriers  in  applied  field  E(t)  and  unperturbed  phonons,  the  second  term  describes  interaction  of 
carriers  with  phonons,  and  the  third  term  describes  carrier-carrier  interaction.  The  Ho  is  given  by 

H0  =1V  (k)CkCk  +Xsh(k)dkdk  +Z^®(q)bqbq 

k  k  k  (17 

+X[MlE(t)e-i‘-'c;d;k  +  M;E(t)ek*'d_lcl  ] 


where  E(t)  =  Eoexp(t2/ip2)  for  a  Gaussian  pulse.  The  electron  and  hole  energies  are  given  by  se(k) 
=  Eg  +  h2k2/2me  and  sh  (k)  =  h2k2/2nih  for  parabolic  bands,  and  by  equation  1  if  the 
nonparabolicity  of  the  conduction  band  is  included.  The  optical  dipole  matrix  element  in  the 
parabolic  bands  approximation  is  given  by 


M 


k 


tr  k2 
2mrEg 


(18) 


which  can  be  approximated  by  a  constant,  Mk  ~  dcv,  for  large  band  gap  values.  The  interband 
electric  dipole  dcv  is  related  to  interband  momentum  matrix  element,  dcv/e  =  hpcv/Egmo,  and  in  the 
effective  mass  approximation  for  two-band  model 


d  /e  = 


-l  +  m0/me 
2Egmo 


-1 1/2 


(19) 


o 

from  which  we  obtain  dcv/e  ~  2.12  A. 


The  carrier-phonon  interaction  Hamiltonian  is  given  by 


Hc-ph  ^[Yk,k+qCk+qkqCk  Tk,k+qCkkq  Ck+q  Yk,k+q^k+qkq^k  Tk,k+q^kkq  ^k+q  ]  ’ 
k,q 


(20) 
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where  for  interaction  with  polar  optical  phonons  the  interaction  matrix  element  yek,k+q  =  -  y\,k+q 
=  MFR(k,k+q)  is  given  in  equation  6,  and 

Hc_c  =  I  V 

k,p,q 


-CkC;cp+qck_q  +  -d(d;dp+qdk_q  -c(d(pd_p+qck_q 


(21) 


where  Vq  is  given  in  equation  8. 

The  kinetic  variables  are  the  distribution  functions  of  carriers  and  phonons  fk“  and  Nq,  and  the 
interband  polarization  pr: 

fl'(t)  =  (c;(t)cl(t))  f;,(t)  =  (d;(t)dl(t))  p1(t)  =  (d.l(t)cl(t))  Nq=(bq*bq),  (22) 


where  all  the  operators  are  in  Heisenberg  representation,  and  the  expectation  values  are  taken 
with  respect  to  initial  state  of  the  system,  which  is  the  vacuum  of  electron-hole  pairs  and  a 
thermal  phonon  distribution.  The  kinetic  equations,  the  SBE,  can  be  derived  using  the  equation 
of  motion  for  a  quantum  operator  Ori(l)  in  Heisenberg  representation,  /ho,6>//  =  [Oh, H]. 


4.1  Single-particle  Problem 

The  evolution  equations  of  non-interacting  carriers  under  the  action  of  Hamiltonian  Ho  describe 
the  generation  of  electron-hole  pairs  in  a  homogeneous  bulk  semiconductor  by  the  action  of  a 
homogeneous  electric  field: 


dt 


df 


dt 


(0) 


=  g£’,(t) 


(23) 


dPk 


i(°) 


dt 


(k)  +  sh( 


k))Pk  +  MkE0(t)e  ‘“pt  (l 


(24) 


with  the  generation  rate 


?<0)(t)  =  ^rMkE(t)e-k°ptP;  -  MkE(t)ek0ptpk 
m  L 


(25) 


These  equations  describe  a  set  of  independent  two-level  systems  driven  by  an  external  field,  with 
a  variable  detuning  Ak  between  the  transition  energy  and  the  light  central  frequency, 

HAk  =8e(k)  +  8h(k)-f7C0p  .  (26) 


The  dynamics  at  this  zero-order  approximation  are  fully  coherent.  The  corresponding  Rabi 
frequency  ( 3 )  is 

Q(R0)E^i  .  (27) 

ft 
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When  the  carrier-carrier  interactions  are  taken  into  account,  this  value  of  the  Rabi  frequency  will 
be  changed  (renormalized)  by  the  Hartree-Fock  terms  (5).  The  intensity  of  the  pulse  is 
proportional  to  |  E(t)  |  2;  therefore,  if  the  width  of  the  intensity  is  tw  =  100  ft,  the  corresponding 
field  pulse  widths  is  xp  =  85  fs.  The  electric  field  E(t)  =  Re  E0exp(/copt  - 1  /xp  ),  and  its  Fourier 
transform  is 


E(M)^e 


-(co-co  )-t  -/4 


(28) 


the  energy  width  is,  therefore,  5E  =  2h/xp  ~  1 5  meV. 


We  have  to  relate  the  field  strength  to  the  power  of  the  incident  pulse  train  used  in  the 
experiments.  The  energy  flux  in  a  single  pulse  is  So(t)  =  SocE2(t),  in  SI  units.  Let  us  assume  that 
the  pulses  come  at  a  repetition  rate  cor.  Then  the  resulting  energy  flux  is 


Sp  =cor  JdtS0(t) 


T- '  ^  I  i  (Q_  T  /2 

®r'Cp£0cEin,0  I  !-e  PP 


(29) 


In  the  experimental  work  described  in  section  6,  the  relevant  wavelength  kp  of  light  is  in  the 
range  370  to  250  nm,  and  the  corresponding  circular  frequency  cop  =  2ttc/a  is  in  the  range  (5.1  - 
7.5)xl015  1/s.  Therefore,  we  can  neglect  the  exponentially  small  term  in  equation  29.  Equating 
the  flux  to  the  power  per  unit  area,  W/As,  where  As  is  approximately  the  area  of  the  sample 
surface  where  the  light  is  absorbed,  we  obtain  a  relation  between  the  electric  field  amplitude  and 
the  power: 


Ein,0  “ 


fn\m 

u. 


W 


s0corXpCAs 


1/2 


(30) 


The  repetition  rate  in  the  experiments  was  250  KHz.  Assume  that  the  incident  spot  area  As  = 
100x100  pm2,  then  for  W  =  100  pW  the  corresponding  value  of  Ein,odcv  is  3.565  meV.  At  normal 
incidence,  the  ratio  of  fields  in  reflected  and  incident  waves  is  (n-l)/(n+l),  and  for  GaN  the 
reflectance  n  =  2.2.  Then  we  obtain  that  in  the  transmitted  wave,  Eo  =  0.927Ein  0  and  Eodcv  = 

3.304  meV,  and  the  corresponding  value  of  the  Rabi  frequency  in  equation  27  is  5.02  THz.  At 
arbitrary  values  of  power  W,  we  obtain  the  following  relation  between  Rabi  frequency  and  the 
power  of  the  incident  pulses  (measured  in  pW): 

Q(r0)  =  0.502  THzx  ^/W/l  pW  .  (31) 

The  equation  24  for  the  polarization  with  the  initial  condition,  pk(-°°)  =  0,  can  be  recast  in  the 
integral  form,  convenient  for  the  discussion  of  the  semiclassical  limit: 

Pk(t)  =  -iGf  e-KOpt  J dxe_iA|cT  E(t  -  x)  [l  -  f ke  (t  -  x)  -  f _hk  (t  -  x)]  .  (32) 

o 
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The  semiclassical  result  is  obtained  by  (1)  adiabatic  approximation  E(t)  — >  E(t)exp(r|t)  with  r\  —> 
0+  and  (2)  Markov  approximation,  in  which  we  assume  that  the  electric  field  and  the  distribution 
functions  in  the  integrand  vary  slowly  compared  to  oscillations  of  the  exponential  factor.  These 
approximations  eliminate  (1)  the  influence  of  initial  correlations  and  (2)  the  dynamic  memory 
effects.  Then  the  polarization  in  equation  32  is  an  instantaneous  function  of  the  carrier 
distribution  and  applied  field.  Using  the  relation 

dxe'(r|+'Ak,T  =  -i  —  +  Jt5(A. ) 

Ak  ’ 

where  P  denotes  the  principal  value  and  5(x)  is  the  Dirac’s  delta-function,  we  obtain  the 
semiclassical  generation  rate  from  equation  25  in  the  form  of  the  Fermi’s  Golden  Rule, 

g">(t)  =  2IIn«2|E(t)|=[l-f1,'(t)-f“k]6(Ak)  .  (33) 


1 


If  the  Markov  approximation  in  equation  32  is  applied  to  the  carrier  distribution  functions  and 
we  retain  the  time  dependence  of  the  Gaussian  field,  we  obtain  the  effects  of  the  non-zero 
spectral  width  of  the  pulse: 


?‘0,(t)=2nf 


Tpe  1  /T|’  [l-fkC(t-T)-f_hk(t~X 


-.00 

)Jjdu 

0 


cos(Akxpu)« 


-(u-t/cp) 


(34) 


2 

For  times  shorter  than  Akxp  ,  the  semiclassical  generation  rate  can  be  approximately  factorized 
into  the  temporal  and  spectral  pulse  shapes  (7): 


r(°>/ 

>k 


(t)  =  2Qf  2xp  [l  -  fke  (t  -  x)  -  f  hk  (t  -  x)]e 


't2  -A?t2/4 

pe  kp  . 


(35) 


4.2  Carrier-phonon  Interaction 

The  contributions  of  carrier-phonon  interaction  are  obtained  by  evaluation  of  the  commutation  of 
carrier  operators  in  equation  22  with  Hc.ph  in  equation  20.  The  resulting  contribution  to  the 
equation  for  the  electron  distribution  function  is  given  by 


dU 


dt 


=  y>  X  fa  (<Vk-q  )  -  Yq  <ck+qbqck  )  +  Y  f  (<b;ck+q 


c-ph 


h  *  / 

-  Yq  (C 


be 

k-quq^k 


(36) 


and  similarly  for  l\  and  pk.  The  lowest  order  approximation,  the  first  order  contribution,  is 
obtained  by  factorizing  the  expectation  values  of  three-operator  products  into  distribution 
functions  and  polarization.  The  resulting  first-order  contributions  vanish  because  (bq  '"j  =  (b+j  =  0 

for  the  eigen-states  of  the  phonon  number  operator,  and,  consequently,  also  in  the  averaging  with 
the  thermal  phonons  density  matrix.  In  order  to  obtain  the  second-order  contribution,  we  have  to 
derive  equations  of  motion  for  the  three-operator  products  on  the  right  hand  side  of  equation  36. 
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Again,  using  commutation  and  anti-commutation  relations  in  the  Heisenberg  equation  of  motion 
for  the  operator  products,  we  obtain 


^  (ckVk-q )  =  [-  se  (k)  +  se  (k  -  q)  +  /to(q)](ckbqc 

+  iz  k  (ckck'ck'+qck_q )  +  yf  (ckdk’dk-+qck_q  ;J 
1/2  k' 


q^k-q 


-  cKb,b,.CM 


})+<(< 


'k^q^q'^k-q+q' 


c+b.bfc,, ,  >  —  ( c  t  ..b'.bc 


'k-q’-q'-q-k-q//j 


(37) 


and  similarly  for  the  other  terms.  The  second-order  approximation  is  obtained  by  factorizing  the 
expectation  values  of  four-operator  products  into  distribution  functions  of  electrons  and  phonons, 
and  polarization,  and  neglecting  the  two-particle  correlations.  Define 


Sq,k,k-q  ~  (CkbqCk-q 


(38) 


and  obtain  the  equation: 


dt 


s;*^,  =  i[-s'(k)  +  e'(k-q)  +  to(q)]s;ki_,  +iTf  IK  +1K‘  -N,fk%]-^r:fk'fk-, 


i  h 

1  e*  * 

-Yq  PkPk-q 


(39) 


\h 


where  we  used  (bqbq )  =  Nq  .  This  equation  has  the  form 


—  S(t)  =  -icoS(t  )  +  Y(t). 
dt 


(40) 


With  the  initial  condition  giving  S(to),  it  can  be  cast  into  the  integral  form 

t-to  (41) 

S(t)  =  S(f0)e-1“(t-t°)  +  J  dxe_1“TY(t  —  t) 

o 

We  perform  the  adiabatic  approximation  by  changing  co  — >  co  +  0+  and  letting  t0  — >  -  oo,  and  also 
eliminate  memory  effects  in  the  scatterings  by  Markov  approximation  in  the  time  integral.  In 
doing  this,  we  assume  that  Y(t)  varies  slowly  as  compared  to  the  exponential  factor  in  the 
integrand.  With  this  we  obtain 

S(t)  =  (£  +  jc8(g>))  Y(t)  =  7iD(co)Y(t),  (42) 

where  T  denotes  principal  value  and 

D(ra)  =  --  +  5(w).  (43) 

ITT  (0 
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The  Markov  approximation  here  relies  on  the  assumption  of  different  time  scales  relevant  for  the 
interaction-induced  correlations  on  one  hand,  and  the  dynamics  of  distribution  functions  and 
polarization  on  the  other  hand,  and  results  in  a  semiclassical  approach  to  the  relaxation  and 
dephasing  processes.  The  limitations  of  this  approximation  for  the  carrier-phonon  kinetics  on  a 
femtosecond  scale  are  discussed  in  reference  1. 


Collecting  all  the  terms  and  using  yeq  =  -  yhq  =  yq  for  interaction  with  polar  phonons,  we  obtain 

(9) 

i(2)  „ 

-  \  fwe'Cphfe  We’Cphfe  1  4-  [pe'Cphn*  pe,cph*  1 

-  “  [Wk— q,lt  k  “  Wk,k— qfk— qj  +  [bk  Pk  ~  Pk  PkJ 


dfk 

dt 


cph 


dfu 

dt 


h,C2) 


cph 


=  -  X  [  <?Ae  -  wkhk^fkh_q]  +  7^[p-£phP-k  -  F-kcph*P-k] 


dpk 


dt 


(2) 


cph 


=  -ZKp-5p--w^Pk_q], 


where  the  internal  polarization  fields  F“’cph  are  given  by 

Fjx.cph  =  iTt^|Yq|2  ^(+pk_q)^(£«  _  £«_q  +  /ia)q), 


(44) 


(45) 


(46) 


(47) 


the  carrier-phonon  scattering  rates  are  given  by 


W, 


a, cph  _  ^TU  \ 


k-q,k 


=  f2> 

+- 


'6(£k  -  £k-q  ±  ^q)  (Nq  +  \  ±  0  -  fk-q)> 


(48) 


and  the  polarization  transition  matrices  are  given  by 


Nq  T 


a=e,h  ±- 


1 

2  +  2)  fk-q 


+  K  +  2±2)6-k%) 


(49) 


The  5-function  part  of  T>  produces  the  scattering  and  dephasing  terms,  while  the  principal  value 
part  produces  energy  renormalization  terms  describing  polaron  corrections  to  the  band  structure. 
It  follows  from  the  last  two  equations  that  in  the  low  excitation  regime,  f“  «  1, 
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(50) 


4.3  Carrier-carrier  Interaction 

The  contributions  of  carrier-carrier  interaction  are  obtained  using  the  Heisenberg  equations  of 
motion.  The  derivation  employs  the  anti-commutation  relations  for  the  electron  and  hole 
operators,  from  which  we  also  derive  the  following  commutation  relations: 

[ck.ck,,ck]  =  5kk'Ck  ,  [ck,ck  ’ckJ  =  — Skkck  ,  [ck.d_k,,ck]=  — 5kk.d_k  ,  [ck,d_k,,ck  J  =  0  . 


The  resulting  c-c  contribution  to  the  equation  for  electron  distribution  function  is  given  by 


iZvj 


dt  \h 

CC  K-tl 


ctctc^c^ )  -  cLcLcwC,  >  -  c'dT  d_yctj, )  +  (ctA\,d_v,c 


k  k  k'+q  k-q  /  X^k-q^k'+q^k'^k  /  \~k~-k’-q~-k'''k-q  /  1  \^k-q“-k'^-k'-q^k 


.  ,  (51) 


and  similarly  for  1* k  and  pk  The  lowest  order  approximation,  the  first  order  contribution,  is 
obtained  by  factorizing  the  expectation  values  of  four-operator  products  into  distribution 
functions  and  polarization.  The  results  are  the  Hartree-Fock  terms: 


FkC(1)Pk  -FkC<1)*Pk 


,  U)  . 

-£■  =Z[siPl+Fir»,(l-fl'-f!'k)-  (53) 

at  in 

CC 

where  Z  is  the  self-energy  due  to  the  e-e  and  h-h  interactions,  and  F  is  the  internal  field  due  to  the 
e-h  interaction  (interband  polarization  field): 

(54) 

k' 

Fr = -2X„.pk. .  (55) 

k' 

The  self-energy  in  equation  54  is  a  real  quantity  that  renormalizes  the  one-particle  carrier 

•  ccf 

energies  ea(k)  and  does  not  contribute  to  a  finite  lifetime.  The  effect  of  the  internal  field  Fk  is 
to  renormalize  the  Rabi  frequency  (3),  and  it  can  be  included  in  the  generation  rate: 

g„W  =gr(t)  +  g(l1|(t)  =  i{[MkE(t)e-1”*'  +F“(1>]P;  -[M’EWe^  +  F“<'>*]pk}  (56) 

The  contribution  due  to  the  internal  field  satisfies  the  sum  rule 


2>i"=o 
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and  does  not  contribute  to  the  generating  of  electron-hole  pairs.  It  changes  the  distribution  of 
carriers  inside  their  respective  bands,  but  in  contrast  to  true  scattering  processes,  it  is  a  coherent 
process  (8).  This  renormalization  adds  nonlinear  terms  to  the  evolution  equations,  which  can  be 
used  as  a  mathematical  model  to  study  optical  nonlinearities  in  semiconductors  (i).  The 
scattering  processes  appear  in  the  second-order  contributions.  These  are  obtained  from  the 
Heisenberg  equations  of  motion  for  the  two-particle  correlations: 

8<CkC£,Ck/+qCk_q>  =  <c£c£,Ck/+qCk_q>  +  f£fk,5k',k-q  (58) 

s  (ckd_k'_qd_k' Ck-q)  —  (^k  d_k'_qd_k' Ck-q)  —  PkPkf ^k',k-q’ 
which  describe  deviations  from  the  corresponding  factorizations. 


The  resulting  equations  for  these  two-particle  correlation  functions  contain  three -particle 
correlations,  which  are  factorized  into  two-particle  correlations,  distribution  functions,  and 
polarizations  (9).  If  the  two-particle  correlations  are  further  factorized  into  distribution  functions, 
the  resulting  equation  for  the  correlation  function  in  equation  58  is  given  by 

ih^L  8<ckck’ck'+qck-q>  =  [— £e(k)  -  Se(k')  +  se(k'  +  q)  +  se(k  -  q)]5(c^cjck'+qck-q)  (59) 

+  (Vq  -  Vk— k"— q)  [fkefke  (l  -  fk\q)  (l  ~  tf-q)  ~  E-qS'+qC1  “  5)tt  ~ 

+  Vq  [pk-qPk  (fk’+q  -  fk  )  +  Pk  +qPk  (fk— q  ~  fk)] 

-  Vk_k'_q  [Pk-qPk'  (fk'+q  -  fk)  +  Pk'+qPk(fk-q  ~  ft)]  ‘ 


This  equation  is  of  the  same  form  as  equation  40  for  carrier-phonon  correlation  function  Se.  We 
transform  it  into  the  integral  equation  as  in  equation  41  and  perform  adiabatic  and  Markov 
approximations.  The  result  is  then  used  in  equation  51  to  obtain  the  second-order  c-c  interaction 
terms  in  the  kinetic  equation  for  electron  distribution  function.  Have  we  retained  the  two-particle 
correlations  in  equation  59,  the  bare  Coulomb  interaction  Vq  =  e‘/(VsoSsq_)  in  the  kinetic 
equation  for  fke  would  become  dynamically  screened  with  the  Lindhard  dielectric  function  of 
electron  gas  ( 1 ,  16): 


(60) 


where 


6r(q,m)  =  1  -  V, 


(o) 


X 


fe 

rk 


fe 

Mc-q 


'k-q 


hoi 


(61) 


In  this  work  we  use  the  static  limit  of  equation  61,  and  for  Vq  in  kinetic  equations,  we  use 
Coulomb  interaction  statically  screened  by  electrons  and  holes,  thus  obtaining  the  screened 
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Coulomb  interaction  in  equation  7  with  the  inverse  screening  length  given  in  equation  8.  The 
dynamic  screening  effects  that  would  lead  to  carrier-plasmon  scattering,  are  excluded  in  this 
approximation.  With  these  approximations,  the  carrier-carrier  interaction  second-order  terms  in 
the  kinetic  equations  are  the  following: 


dt 


(2) 


=  -XK'V*  -  WS-qfk-q]  +  E  KCCPi-  “  TC‘pJ 


dt 


cc 


w. 


h.cc 


fe 

q.^k 


w, 


h.cc 


k,k-q 


h,cc(2) 

-k 


P-k 


dPk 


dt 


(2) 


=  -BwkP-q.kPk-w£k'!qPk_q], 


„h,cc(2)* 

F-k 


(62) 


(63) 


(64) 


where  the  internal  polarization  fields  F"‘ct(2’  are  given  by 


i-.a,cc(2) 

Fk 


4-  oa 

+  £k' 


'k'+q 


)■ 


(65) 


the  carrier-carrier  transition  matrices  are  given  by 

\  |V,|2  X  Z  ®6“k-q  -  4  +  <£+,  -  4)(i  -  «-„)  [«'  (l  -  fk"'+q) 

at  kf 

-Pk'+qPkf]+  C.C., 


wkai“,k 


(66) 


and  the  polarization  transition  matrices  are  given  by 


Wk-q,k  =  \  lvq|2  X  Z  -  £k  +  £k'+q  -  £k0  [-Pk'+qPk' 

a, at  k' 

+  fk/  (l  -  fCq)  (1  -  fk— q)  +  fk“-qfk“:+q(l  ~  «0]  • 


(67) 


4.4  Kinetic  Equations  to  the  Second-order  in  Carrier-phonon  and  Carrier-carrier 
Interactions 

With  the  interaction  terms  derived  in  section  4.3,  the  equations  of  motion  can  be  written  as 

ira 

=  gkW  -  £[wk“.q.kfk“  -  wk“k_qfk“_q] 

q 
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(69) 


^  ^  (4e  +  £t k)pk  +  2 rkp(t)(l  -  fk»  -  fkh)  -  Y  KP-q,kPk  -  wpk_,Pk_q], 

q 

with  the  renormalized  energies  of  carriers  (Hartree-Fock  terms) 

££  =  el  +  Zk,  (70) 

the  fields  given  by  the  sum  of  the  applied  field  the  and  internal  polarization  fields 

=  MkECOe-10^1  +  Fk  cph  +  f“’cc(1)  +  f“’cc(2) 

T?  =  +  F“’cc(1), 

the  generation  rate 

(73) 

and  with  transition  matrices  Wkk_q  and  Wkk  being  the  sums  of  the  second-order  carrier- 
phonon  and  carrier-carrier  interaction  transition  matrices.  The  transition  matrices  Wkk_q  (a  =  e, 
h)  for  the  distribution  functions  are  real  quantities  but  are  not  necessarily  positive.  The  transition 
matrix  for  the  polarization  Wk  k  is  a  complex  quantity.  Its  real  part  describes  scattering 

processes  that  contribute  to  the  dephasing  of  polarization,  and  the  imaginary  part  describes 
second-order  contributions  to  the  band  gap  renormalization.  The  second-order  energy  shifts  are 
usually  negligible  as  compared  to  the  Hartree-Fock  energy  renormalization  ( 1 ). 

4.5  Luminescence  Spectrum 

In  this  work,  we  will  not  consider  the  emission  of  coherent  photons  possible  in  the  presence  of 
coherent  interband  polarization  (8),  and  the  luminescence  spectrum  will  be  evaluated  from  the 
Fermi  Golden  Rule.  It  results  from  evaluation  of  the  rate  of  emitted  photons  using  Markov 
approximation  of  the  carrier-photon  interaction  in  the  luminescence  to  the  second-order.  With 
the  carrier-carrier  distribution  functions  obtained  from  the  kinetic  equations  68  and  69,  and  the 
momentum  of  emitted  photons  approximated  by  zero,  the  luminescence  spectrum  is  given  by 

lem(w)  =  +  £"k  ~  n(°)fkfk  '  ^ 

k 

Thus,  the  luminescence  spectrum  is  proportional  to  the  product  of  electron  and  hole  distribution 
functions  multiplied  by  the  joint  density  of  states.  The  numerical  solutions  obtained  using  a 
generalized  Monte  Carlo  simulation  (8,  9)  will  be  shown  in  section  7. 


(71) 

(72) 


18 


5.  Relaxation  Time  Approximation 


A  relaxation  time  approximation  of  collision  terms  in  kinetic  equations  allows  use  of  a  simple 
model  to  study  coherent  effects  related  to  optical  interband  transitions.  In  this  model,  different 
momentum  values  decouple,  and  one  ends  with  effectively  a  two-level  system  whose  dynamics 
can  be  described  by  evolution  of  a  three  component  Bloch  vector  (5).  Define  the  components  as 
follows: 

-ico  t  sb  ico  t 

Ulk(t)  =  e  Ppk(t)  +  e  Ppk(t) 


(  -ico  t  sk  ico  t 

U2k(t)  =  -i  e  Ppk(t)-e  Ppk(t) 


U3k  =  1  -  f k  (t)  -  f k  (t) 


and  define  detuning  by  equation  26:  /iAk  =  £k  +  £k  —  /to)p  .  Then  from  equations  23-25,  and 
allowing  for  relaxation,  the  optical  Bloch  equations  are 


^  =  _AU 

dt  k  2k  T2k 


^  =  AkUlk+2<E(t)U3k--^ 


=  -2D^E(t)U2k - ^ 


with  the  initial  conditions 


U3k(-0))=1,  Ulk(— oo)  —  U2k(— oo)  =  0 


The  equations  without  relaxation  describe  a  rotation  of  vector  Uk  about  vector  =  [-2f2R(0,E(t), 
0,  Ak],  in  which  case  they  can  be  written  as  one  vector  equation:  dU/dt  =  H  x  U  .  In  view  of  the 
discussion  in  section  3,  we  assume  that  the  relaxation  in  this  model  is  due  to  carrier-phonon 
interaction,  and  we  set  the  population  decay  time  T \  (“longitudinal  relaxation  time”)  to  be  given 
by  the  inverse  of  the  sum  of  electron-phonon  and  hole-phonon  rates: 


pe-ph  .  ph-ph  1 


In  view  of  the  approximation  in  equation  50,  the  polarization  decay  time  T2  (“transverse 
relaxation  time”)  is  twice  the  population  decay  time,  T2k  =  2Tik. 
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In  the  absence  of  relaxation,  the  carrier  distribution  function  obtained  from  the  optical  Bloch 
equations  will  exhibit  time  oscillations  with  the  Rabi  frequency.  If  the  difference  between  the 
pulse  central  frequency  and  the  band  gap  is  appreciatively  smaller  than  the  LO  polar  phonon 
frequency,  one  expects  a  significantly  reduced  carrier-phonon  relaxation  in  GaN  at  room 
temperature,  in  which  case  the  Rabi  oscillations  should  be  seen  in  the  time-resolved 
luminescence.  The  luminescence  evaluated  in  this  model  at  room  temperature  is  shown  in  figure 
7  for  three  values  of  excitation  energy  relative  to  the  band  gap,  at  two  different  values  of 
excitation  power. 


Time  (ps) 


Figure  7.  Luminescence  evaluated  in  relaxation  model  at  room  temperature  for  three  values  of 
excitation  energy  relative  to  the  band  gap,  at  two  different  values  of  excitation  power. 

These  results  indicate  a  possibility  of  observing  coherent  oscillations  in  the  time -resolved 
luminescence  experiments  in  high  quality  GaN. 


6.  Experimental  Studies  of  the  Time  Resolved  Spectroscopy  of  UV 
Luminescence 


A  schematic  representation  of  the  experimental  setup  (7)  is  shown  in  figure  8.  A  250-KHz 
regeneratively  amplified  Ti:Sappfire  laser  produces  150  fs  pulses  at  800  nm  wavelength.  Part  of 
this  beam  is  split  off,  acting  as  a  gating  pulse,  and  the  other  part  pumps  an  optical  parametric 
amplifier  tunable  from  720  to  450  nm.  Output  from  the  amplifier  is  again  doubled  in  a  nonlinear 
optical  crystal  (2Xtal)  to  produce  pulses  tunable  from  360  to  225  nm.  The  pump  pulse, 
compressed  to  a  width  of  approximately  80  fs,  is  focused  onto  the  sample  where  it  photo-excites 
electron-hole  pairs.  The  luminescence  from  the  sample  is  collected  using  non-dispersive  optics 
and  focused  on  another  nonlinear  crystal  (NLXtal),  where  it  is  mixed  with  the  gating  pulse  to 
yield  a  down-converted  signal  at  lower  frequency  (visible  light)  that  is  spatially  separated  from 
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the  other  beam.  The  down-converted  signal  is  passed  through  a  spectrometer  and  detected  with  a 
photomultiplier  tube  in  photon  counting  mode.  The  sample  used  in  these  experiments  was  a  1  pm 
thick  homoepitaxial  GaN  film  deposited  on  a  GaN  substrate.  The  arrival  time  of  the  gate  pulse 
relative  to  the  photoluminescence  is  used  to  follow  the  temporal  evolution  of  UV 
photoluminescence  with  350  fs  resolution. 


Figure  8.  A  schematic  representation  of  the  setup  for  the  time-resolved 
photoluminescence  experiment.  A  Ti:Sapphire  laser  produces 
150  fs,  5  mJ,  800  nm  pulses.  Part  of  the  beam  pumps  an  optical 
parametric  amplifier  tunable  from  720  to  450  nm.  An  up-converted 
beam  with  pulses  tunable  from  360  to  225  nm  and  widths  of  50  to  80  fs, 
pumps  energy  from  3.44  to  5.5  eV,  suitable  for  studying  nitride  materials, 
generates  electron-hole  pairs  in  a  GaN  sample. 

The  experiments  were  performed  at  two  values  of  pump  pulse  power,  0.66  mW  and  2  mW, 
generating  about  8x  1018  and  2.5  x  1019  electron-hole  pairs,  respectively.  In  the  experiments, 
the  luminescence  was  measured  at  a  fixed  energy,  thus  allowing  us  to  watch  the  carrier 
distribution  pass  through  a  spectral  window,  as  shown  in  figure  9. 
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Figure  9.  The  upper  part  illustrates  the  energy  position  of  2  mW  pump  pulse,  at  308  nm 

wavelength,  and  four  different  values  of  emission  energy.  The  band  edge  exciton 
luminescence  lineshape  is  also  shown.  The  lower  part  shows  the  time 
dependence  of  the  measured  photoluminescence  at  those  values  of  emission 
energy,  set  near  the  band  edge  (black),  hmLo  (red),  2  ho)|  0  (green),  and  3htt>LO 
(blue),  where  ha>LO  referes  to  the  energy  of  a  polar  longitudinal  phonon,  91.2  meV. 

Room  temperature  time-resolved  photoluminescence  obtained  with  a  308  nm  pump  (about 
75  meV  above  the  conduction  band  edge)  is  also  shown  in  figure  10  for  two  values  of  the 
excitation  power,  0.66  mW  and  2  mW. 


Figure  10.  Time  evolution  of  photo-luminescence  from  GaN  at  different  wavelengths.  The  pump  wavelength 
is  308  nm,  at  two  values  of  the  pump  power,  0.66  mW  (left  panel)  and  2  mW  (right  panel). 
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The  data  are  indicative  of  the  rapid  formation  of  a  non-thermal  distribution  associated  with  the 
interplay  of  carrier-carrier  scattering  (electron-electron  and  hole-hole)  and  electron-LO  phonon 
scattering  during  carrier  injection.  After  injection,  a  measurable  rise  time  to  peak  emission  is 
present  and  is  a  function  of  excitation  energy,  carrier  density,  and  emission  energy.  Since  the 
time-resolved  signal  at  a  given  wavelength  is  proportional  to  the  product  of  the  electron  and  hole 
distribution  functions  at  a  given  momentum,  these  plots  provide  a  measure  of  the  evolution  of  the 
carrier  distributions  as  they  cool  down  to  the  band  edge  and  into  possible  deeper  band  tail  states. 

Carrier  density-dependent  data,  shown  in  figure  11,  shows  the  effects  of  screening  in  evolution 
of  electron-hole  distribution.  For  an  initial  excess  energy  of  0.75  eV  (308  nm  excitation  pulse),  a 
shift  of  the  355.6  nm  (1LO  phonon  energy  above  the  conduction  band  edge)  emission  peak  from 
0.9  ps  to  2.2  ps  is  seen  for  a  change  in  density  from  8x10  to  2.5x10  cm  "  (660  pW  to  2  mW 
pump).  The  fact  that  the  emission  peak  is  reached  later  at  this  wavelength  at  higher  carrier 
density  implies  that  the  electrons  relax  more  slowly  due  to  more  efficient  screening  of  the 
electron-LO  phonon  interaction. 


Figure  11.  Pump  power  dependence  of  355.6  nm  emission 
showing  an  increase  in  rise  time  at  higher  density 
due  to  screening  of  the  electron-LO  phonon 
interaction. 

These  experiments  allowed  us  to  study  initial  thermalization  and  hot  carrier  dynamics  in  photo- 
excited  GaN.  Screening  effects  were  observed  as  a  function  of  excitation  intensity,  and  a 
dependence  on  the  excitation  photon  energy  is  seen  in  the  dynamics  of  the  photo-excited  carrier 
distributions. 
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7.  Monte  Carlo  Simulation  of  Kinetic  Equations 


The  numerical  method  we  used  in  the  simulation  is  the  Monte  Carlo  method  generalized  by 
Kuhn  and  Rossi  to  account  for  the  effects  of  the  interband  polarization  (8).  The  various 
contributions  to  the  carrier  dynamics  in  equations  68  and  69  can  be  separated  into  two  parts,  the 
coherent  part  and  dephasing  part.  Let  “Rvk  be  a  kinetic  variable,  being  fke,h  if  v  =  e,h  or  pk  if  v  =  p. 
The  kinetic  equations  can  be  schematically  written  as 


—  7?v  -  —  7?v 
dt  k  "  df  k 


coher  UL 


incoh 


,  v  =  e,  h,  p 


(79) 


with 


coher 


Sk°({32v})  +  ^5k)(WV}) 

1*0 


incoh 


j*' 


(80) 


where  is  a  functional  of  kinetic  variables,  and  is  the  scattering  rate  associated  with 
kinetic  variable  !RV  for  a  transition  k  — »  k'  induced  by  the  jth  interaction  mechanism.  These  rates 
are  real  positive  quantities,  unlike  the  full  transition  matrices  W^,. 

The  coherent  part  describes  the  generation  of  carriers  and  polarization,  which  are  highly 
correlated  due  to  the  external  field.  This  part  is  evaluated  by  a  direct  numerical  integration. 
Under  the  approximations  used  in  the  derivation  of  kinetic  equations  in  this  work,  the  incoherent 
part  has  basically  the  same  form  as  used  in  traditional  Monte  Carlo  simulation.  The  only 
difference  is  that  the  polarization  is  a  complex  quantity. 

The  total  time  is  divided  into  time  steps.  At,  and  the  system  is  initialized  in  vacuum  state  at  to.  In 
this  state,  the  distribution  functions  and  the  polarization  are  zero.  For  each  step,  we  compute  the 
coherent  part  by  the  direct  integration.  In  the  Monte  Carlo  “sampling”  of  incoherent 
contributions  for  each  kinetic  variable  !RV,  we  use  an  ensemble  of  N  particles  generated 
randomly  according  to  the  distribution  function  1 77y | .  For  the  polarization,  a  phase  is  attached 
according  to  the  phase  of  pk.  The  traditional  Monte  Carlo  simulation  is  performed  according  to 
the  scattering  rates  W^,.  A  more  detailed  description  can  be  found  in  reference  9.  The 
conduction  band  nonparabolicity  and  static  screening  were  included  in  the  calculations.  Once  the 
distribution  functions  are  obtained,  the  luminescence  is  evaluated  from  equation  74.  The  results 
in  figure  12  show  that  the  rise  of  the  experimentally  observed  luminescence  and  the  associated 
rise-times  are  in  reasonable  agreement  with  experimental  results  in  section  6,  confirming  that  the 
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carrier-polar  optical  phonon  interaction  is  the  dominant  scattering  mechanism  in  the 
subpicosecond  evolution  of  photo-excited  carriers. 


Figure  12.  Theoretical  luminescence  from  photo-excited  GaN  obtained  by  Monte  Carlo  simulation  of  kinetics 
equations.  Luminescence  is  shown  as  function  of  time  at  four  different  energies,  at  the  band  edge 
and  at  one,  two,  and  three  times  LO-phonon  energy  above  the  band  gap. 


8.  Conclusions 


We  have  applied  the  SBE  formalism  to  study  the  subpicosecond  dynamics  of  photo-excited 
carriers  in  GaN.  The  kinetics  equations  were  derived  using  a  two-band  model  and  taking  into 
account  the  coherent  interaction  of  semiconductor  with  the  short  laser  pulse.  Various  scattering 
terms  due  to  interaction  of  carriers  with  polar  optical  phonons  and  carrier-carrier  interaction  were 
taken  into  account.  Effects  of  the  conduction  band  nonparabolicity  and  static  screening  were 
included  in  our  theory.  Prior  to  the  evaluation  of  kinetics,  we  calculated  scattering  rates  and 
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found  that  the  relaxation  of  carriers  in  bulk  GaN  is  dominated  by  electron-polar  phonon 
scattering,  with  rates  much  higher  than  the  corresponding  rates  in  GaAs,  for  example.  We  also 
found  that  the  plasmon  emission  by  energetic  electrons  can  be  a  substantial  contributing 
mechanism  to  the  hot  electron  relaxation  on  the  subpicosecond  scale,  though  not  nearly  as 
important  as  carrier-optical  phonon  scattering.  Further  investigation  of  possible  plasmon  effects 
may  be  included  in  the  future  theoretical  work  on  hot  carrier  dynamics  in  GaN.  This  will  involve 
using  dynamic  carrier  screening  rather  than  the  static  approximation.  We  find  reasonable 
agreement  with  the  experimental  results  on  the  time-resolved  photoluminescence  rise-time 
dependence  on  excitation  energy  and  carrier  density,  thus  confirming  the  dominant  role  of  polar 
optical  phonons  in  the  subpicosecond  dynamics,  and  the  effects  of  screening.  The  results  of  the 
relaxation  model  suggest  a  possibility  of  observing  coherent  polarization  effects,  the  Rabi 
oscillations  in  luminescence,  possibly  using  different  wide-band-gap  polar  material  samples  at 
lower  pump  power.  An  extension  of  the  experiments  to  lower  temperatures  would  also  be  of 
interest  to  further  clarify  the  role  of  polar  phonons  in  subpicosecond  hot  electron  dynamics. 
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List  of  Acronyms  and  Abbreviations 


2Xtal 

nonlinear  optical  crystal 

AlGaN 

aluminum  gallium  nitride 

Fs 

femtosecond 

GaAs 

gallium  arsenide 

GaN 

gallium  nitride 

LO 

longitudinal  optical 

NLXtal 

nonlinear  crystal 

RPA 

Random  Phase  Approximation 

SBE 

Semiconductor  Bloch  Equation 

TRPL 

time-resolved  photoluminescence 

UV 

ultraviolet 
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